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Abstract 

Hadron bubbles that nucleate with radius R nuc in a quark sea (if the phase transition 
is first order) are shown to be unstable to the growth of nonspherical structure when the 
bubble radii exceed a critical size of 20 — 10 3 R n uc- This instability is driven by a very thin 
layer of slowly diffusing excess baryon number that forms on the surface of the bubble 
wall and limits its growth. This instability occurs on a shorter length scale than those 
studied previously and these effects can thus be important for both cosmology and heavy 
ion collisions. 
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If a first order phase transition occurs when quarks are confined into hadrons at a 
temperature Tq = 100 - 200 MeV, many interesting physical consequences can arise [1]. 
Previous investigations [1, 2] have shown that the baryon number density can be higher 
in the quark phase than in the hadron phase. In the case of cosmology, this asymmetry 
can lead to inhomogeneous distributions of baryons after the phase transition [1, 3]. This 
possibility has led to numerous studies of big bang nucleosynthesis in the presence of 
baryon inhomogeneities [3 - 7] . Another possible relic of such a phase transition is stable 
macroscopic quark matter - quark nuggets - which may survive to the present epoch 
and contribute to the dark matter of the universe [1, 8]. The effects of a first order 
quark/hadron (Q/H) phase transition are also, in principle, detectable in relativistic 
heavy ion collisions. 

In this Letter, we assume that the Q/H phase transition is first order and study a 
new instability in the growth of hadron bubbles. When the background plasma, either 
in the early universe or in a heavy ion collision, cools to a temperature below the critical 
temperature Tc, bubbles of the hadronic phase nucleate and begin growing spherically. 
Witten [1] has argued that baryon number prefers to be in the quark phase. As a result, 
the bubble rejects baryon number as it grows. Because baryon number diffuses very 
slowly, it is confined and concentrated into a thin layer on the surface of the bubble. 
We argue that this boundary layer must limit the growth of the bubble and cause it to 
become morphologically unstable. 

In this analysis we assume that the latent heat of the phase transition is transported 
radiatively, mostly by neutrinos, since they have the longest mean free path. The rela- 
tive importance of radiation and hydrodynamic flow in removing the latent heat remains 
unsettled [1, 4, 9, 10]. Previous work has studied the (stable) evolution of hadron bub- 
bles when hydrodynamic flow dominates the removal of latent heat [11, 12]. However, 
radiative processes are present and perhaps dominant from the moment of bubble nucle- 
ation onward. In addition, we expect that the new instabilities (due to baryon number 
diffusion) presented here can arise regardless of the type of heat flow. 

In previous work [9] , we showed that when radiative effects dominate, hadron bubbles 
are unstable for bubble radii greater than the mean free path of the neutrino (~10 cm) 
because the heat flow becomes diffusive at that size scale. Hydrodynamic instabilities 
during hadron bubble growth [13], if present, occur on a smaller scale of order 10~ 6 
cm. The instability proposed here occurs on yet smaller scales, at most 10 5 fm, and 
thus greatly alters the manner in which the phase transition proceeds. The resulting 
inhomogeneities in baryon number are too small (in length scale) to affect big bang 
nucleosynthesis. In fact, this size scale is small enough that this new instability can 
occur in heavy ion collisions, whereas the previously studied instabilities cannot. We 
emphasize, however, that significant uncertainties exist in both the QCD physics and our 
knowledge of nonequilibrium pattern formation. 

We start by considering a locally flat section of a bubble wall at which the hadronic 
phase is advancing at speed v into the quark sea. Previous work [1, 2] has shown that, 
when these phases are in equilibrium with each other, the baryon number in the quark 
phase n q can be much higher than that in the hadron phase n^; the ratio rih/n q is 
typically 0.003 - 0.2. For simplicity, we assume that this ratio vanishes and henceforth 
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consider only quantities in the quark phase. 

In a frame of reference moving with the bubble wall, the diffusion equation for the 
baryon number density n in the quark phase is 

where the direction of motion is along the z-axis and where D is the diffusion constant. 
We will assume that v ~ constant; we justify this assumption below. Since quarks carry 
the baryon number in this phase, the diffusion constant is given by D = (v q )£ q /3, where 
(v q ) is the mean speed of the quarks ((v q ) is fairly close to the speed of light c) and £ q is 
the mean free path of a quark (about 1 fm for Tc ~ 200 MeV) . 

Deep in the quark phase, very far from the advancing hadronic front, the baryon 
number density has a value appropriate for thermal equilibrium at temperature T^, i.e., 
^oo = g (fi/T) oo T^, where the constant g depends on the number of degrees of freedom 
in the quark phase and where fi is the chemical potential. For our universe, standard big 
bang nucleosynthesis [14, 15] constrains \i to be very small, i.e., (ji/T)^ ~ 10~ 9 — 10~ 8 . 
For heavy ion collisions, (n/T) 10" 4 - 10" 2 . 

The diffusion equation (1) must be supplemented by boundary conditions at the 
Q/H interface. We assume local thermodynamic equilibrium, so that the baryon num- 
ber density nj and temperature Tj at the interface are related by a Gibbs-Thompson 
condition [16, 17] which, in this case, takes the form 

T I = T c (l-d K)~/3T A nj. (2) 

Here, k is the interfacial curvature; do = aTcC p /L 2 is a capillary length and is propor- 
tional to the surface tension cr; the specific heat is C p = 0(Tq); and the latent heat is 
L = 0(T C ). Lattice gauge theory calculations [18] have estimated the surface tension to 
be a ~ 0.1 Tq. The third term in Eq. (2), f3T A n 2 j, is the decrease in coexistence temper- 
ature due to the presence of baryon number (solute). Notice that this term is quadratic 
rather than linear in the solute concentration nj. We have defined a convenient energy 
scale Ta = L/C p , which we expect to be 0{Tc). To estimate the value of /?, we use 
a bag model for the equation of state (and we take nh/n q = 0). We find that (3 ~ 0.6 
(Tc /Ta)T^ 6 ~ 0.6 TJ 6 , where we have used the results of Ref. [11] in obtaining the 
numerical coefficient. 



Continuity of baryon number at the interface requires 



njv i = —D 



v ■ Vn 



2 = 2/ 



(3) 



where zj is the interface position, v is the unit normal to the interface, and vj is its speed 
in the z/-direction. The physical meaning of Eq. (3) is that the rate at which excess 
baryon number is rejected from the interior of the hadronic phase is balanced by the rate 
at which baryon number diffuses ahead of the two-phase interface. This effect produces 
the destabilizing boundary layer of excess baryon number on the surface of the growing 
bubble. 
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To complete the statement of our model, we must impose a condition of heat continu- 
ity at the moving interface. We can consider a full treatment of the heat transport either 
through a diffusion equation for temperature or through a radiative transfer equation [9] . 
For the parameter regime of interest, however, neutrinos carry away most of the energy 
and the mean free path of the neutrino is ~ 10 cm, much greater than the size scale of 
the bubble. Hence, the condition of thermal continuity reduces to the simple form [9] 

vo = c ^fr — • ( 4 ) 

Our assumption here is that all of the latent heat produced in the phase transition is 
efficiently carried out of the system by radiation and that the interface moves at the fastest 
speed for which this condition holds. [We have implicitly assumed that the bubbles are 
far enough apart that the latent heat released during the evolution of a given bubble 
does not affect other bubbles. If, however, the mean bubble separation is less than the 
neutrino mean free path, mutual heating effects should be considered.] 

We now consider the stability of the moving bubble wall. As we will see, the situ- 
ations of interest to us are those for which ni ^> rioo. In this case, the baryon number 
density in the boundary layer ahead of the unperturbed wall at zj = is accurately given 
by Eqs. (1) and (3) to be ni exp(— 2z/£), where the layer thickness £ = 2D/v is very 
small for small D. To study the stability of this wall against periodic deformations of 
wave-vector k, we consider linear perturbations of the form 

Sn(z, x, t) ~ h k exp [zk • x — q k z + u k t\ , 

6zi(x., t) ~ Zk exp [zk • x + u k t\ , (5) 
5T/(x, t) ~ T k exp [zk ■ x + uj k t\ , 

where x denotes the position in the plane of the wall. After performing a standard 
stability analysis [19], we find that the amplification rate uj k is 



v(q k £-2)/£ 



1 + D(q k £- 2) /2(3njc£ 



_ d£k< 



A(3nj 



(6) 



where d = (T c /T A )d and q k £ = 1 +(1 + k 2 £ 2 + uj k £ 2 /D) 1 / 2 . The crucial result is that 
the wall is unstable, i.e., u k > 0, for wave- vectors k < he, where 

k c = {4/3nj/d£) 1/2 . (7) 



We now propose two separate estimates of the characteristic length scales for patterns 
produced by this instability. First, consider a growing hadronic bubble with radius R 
and v = R. If we assume that all the baryon number that was swept out by the growing 
bubble remains in the boundary layer at its surface, then 2nR 2 £ni = 4nR s n 00 /3. With 
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this assumption, Eqs. (2) and (4) can be combined to yield a differential equation for 

v R . 2d Q 2 (RR\* 

c=- = A -R-^{3d) ' (8) 
where A = {Tq — T' 00 )/Ta is a dimensionless measure of the undercooling. The condition 
R — determines the nucleation radius R nuc = 2d/ A. For Tc = 200 MeV and a = 0.1T C , 
Rnuc ~ (0.2 fm)/A. The relevant value of A is unknown, although cosmological studies 
usually assume small values, e.g., A ~ 10 -2 - 10~ 4 ; for heavy ion collisions, we assume 
a larger value A ~ 1. The bubble becomes unstable for R 3> R nuc - When the instability 
occurs, however, the third term in Eq. (8) is still very small (the results presented below 
verify this claim); thus, our implicit assumption of a nearly constant growth speed v is 
valid. 

Our stability analysis implies that the bubble becomes morphologically unstable 
when its radius is appreciably larger than Rc = 2n/kc- (This result can easily be 
confirmed by a more systematic analysis in spherical coordinates.) Using our formula for 
kc and the conservation of baryon number assumption to eliminate nj, we find 
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(9) 



Notice that D/cd is of order unity. For cosmology, /Sn 2 ^ ~ 10 -17 . If we take an un- 
dercooling of A = 10 -3 , we find Rc/Rnuc ~ 10 3 , or Rc ~ 10 5 fm. [Notice that our 
approximation ni/noo ^ 1 is valid; for this case, ni/noo = 2R/3£ ~ 10 3 .] For heavy ion 
collisions, we obtain Rc / R nuc ~ 

20, or R c ~ 4 fm. 

The problem with the above argument is that it says little about the final con- 
figuration of the bubble. Indeed, linear stability analysis is notoriously unreliable for 
predicting quantitatively the final relevant length scales for such highly nonlinear pro- 
cesses. Unfortunately, the present situation is not, as far as we know, directly analogous 
to pattern-forming systems that have been studied experimentally. The closest analogy 
that comes to mind is directional solidification of a dilute alloy where a planar front be- 
comes unstable and breaks up into elongated cellular structures. The excess solute (here, 
baryon number) that accumulates ahead of the front diffuses along the boundary layer 
and ultimately is trapped, in concentrated form, in the narrow, sometimes filamentary, 
interstices between the cells. The filaments themselves are usually unstable and break up 
into chains of droplets, which are made up of concentrated solute and trail behind the 
moving front. Although this mechanism is well known in metallurgy, no reliable theory 
for predicting the cellular spacing is available. 

As our second estimate of characteristic length scales produced by the instability, 
we propose the following. Suppose that the Q/H interface leaves behind it baryon-rich 
filaments with characteristic spacing W. Clearly, W cannot be larger than the stability 
length 2n/kc] we therefore assume marginal stability and set W ~ k^, 1 . Let us further 
assume that the radii of the filamentary interstices, where they join the moving interface, 
are roughly equal to the diffusion length t or, equivalently, i? nMC , both of which are 
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about the same size and seem to be the only relevant length scales. (We expect that this 
assumption is the weakest part of our argument. In any case, further work should be 
done to check its validity.) We can now use global conservation of baryon number, i.e., 
riooW 2 ~ nit 2 . From this result and our formula for kc, we find 



W 



R r 



pn 



2 

oo 



1/6 

(10) 



which implies cosmological values of W about one order of magnitude smaller than those 
predicted by Eq. (9). 

We can take this argument one step further and estimate the size of the filaments far 
behind the advancing front. The filaments will contract to a radius w where the baryon 
number density n/ is sufficiently large that T = and v = 0; the Gibbs Thompson 
relation applied on a tubular surface is 

A + rf/ro = pn 2 f , (11) 

where the curvature has changed sign (from Eq. [2]) because the hadronic phase is now 
on the outside. Conservation of baryon number requires that nil 2 ~ n/ro 2 , and we thus 
obtain 

1/9 

(12) 



VJ 



R 



nuc 



pn 



which implies that w is of order 1 fm. Thus, the size scale of any baryon number 
inhomogeneities produced in the phase transition will be much smaller than that required 
(> 100 cm) to affect big bang nucleosynthesis. 

In summary, we have shown that the growth of hadron bubbles during the Q/H 
phase transition can be unstable due to a new instability mechanism. The instability 
first sets in when the bubble grows to a radius of ~ 20 — 10 3 times the original nucleation 
radius. These instabilities can greatly alter the long term evolution of the bubbles. We 
have presented here a marginal stability hypothesis in which the moving interface leaves 
behind tubular regions which are rich in baryon number and which ultimately break up 
into droplets of size w ~ 1 fm. 

The most important of our assumptions are as follows: (A) We have assumed that 
the Q/H phase transition is in fact first order. We stress that the evidence on this 
issue remains divided. (B) We have assumed local thermodynamic equilibrium and have 
written the thermodynamic condition at the bubble interface as the Gibbs Thompson 
relation given in Eq. (3). The true form of this relation depends on the details of 
QCD physics which are not yet well understood. (C) We have assumed that hadron 
bubble growth is limited by the diffusion of baryon number away from the interface. 
This assumption is expected to be valid for sufficiently large baryon number densities 
(or, equivalently, large u/T) and/or sufficiently large bubble radii. We note that the first 
of these conditions is more likely to be met in heavy ion collisions, whereas the second is 
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more likely to be met in the early universe. (D) We have assumed that all of the latent heat 
produced by the phase trasition is carried away by radiation and hence no hydrodynamic 
flow occurs. We note that even when this assumption is violated, the requirement of 
baryon number diffusion can still lead to instabilities which are qualitatively similar to 
those presented in this paper. (E) We have invoked a marginal stability hypothesis 
to describe the long term evolution of the system. This type of behavior is not well 
understood even in laboratory systems, much less in the early universe. 
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